Exact solution for the two-site Kondo-Lattice Model 
- a limiting case for an insulator 
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The Kondo-lattice model is well established as a method to describe an exchange coupling be- 
tween single conduction electrons and localized magnetic moments. As a nontrivial exact result the 
zero-bandwidth limit [atomic limit) can be used to test approximations for this model. As soon 
as the translational symmetry is broken (for instance by sublattice structures) it is necessary to 
consider more than one lattice site. Therefore, we study as a starting point for generalizations the 
situation of a two-site cluster. An equation-of-motion approach is chosen to obtain the one-particle 
Green's function. In order to determine the spectral weights of its energy poles, we derive different 
possibilities for the calculation of the involved correlation functions. In this paper the analytical 
exact result for the situation of an insulator is presented. In a forthcoming article we generalize the 
calculations to arbitrary electron densities, keeping as the only constraint S = 1/2. 
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I. INTRODUCTION 



The. idea of the Kondo lattice model goes back to 
Zeneiu, one of the first scientists who tried to give a 
qualitative explanation of the physical behaviour of tran- 
sition metals. His model is based on the existence of a 
local magnetic moment Sj per lattice site i composed 
by electrons of partially filled electron shells. This as- 
sumption is well fulfilled in the 3d and 4f shells of tran- 
sition metals and lanthanides respectively, since these 
shells are screened very effectively by outer filled electron 
shells. A strong Hund's rule coupling within the shells 
in connection with the notion of the Wigner-Eckart the- 
orem leads to a well defined quantum number S given 
by Sf = h 2 S(S + 1). Due to exchange mechanisms there 
exists an effective coupling between the localized spins, 
which has the Heisenberg form 



(1) 



The electronic behaviour of these materials is de- 
termined by another class of electrons, namely the 5s 
and 5d/6s shells of the 3d- and 4f-systems, respectively. 
These electrons are itinerant and their uncorrelated prop- 
agation through the lattice is described by the Hamilto- 
nian 



(2) 



Using second quantization, the Fermi operator c J(T an- 
nihilates an electron with spin a at site j whereas c\ a 
creates one at site i. The hopping integrals X^- are con- 
nected by Fourier transformation to the single electron 
Bloch energies 



^ = ^E< k ) eik(Ri ~ R3) - 



(3) 



Zener pointed out that for the transition metals it is 
insufficient to handle the magnetic and the electronic par- 
tial system independently. Instead, he suggested an on- 
site interaction of the itinerant electron spin er with the 
localized magnetic moments Si of the form 



K 



(4) 



Due to the close relation to Kondo systems, the Hamil- 
tonian Ti. s + Ti s f is often called Kondo lattice mode/a. 
Another common name is s-f model or s-d model em- 
phasizing the class of described materials. In the case 
of manganites the same Hamiltonian is discussed in the 
limit Jh 3> Tij and called double exchange model, based 
on a special mechanism for the hopping of electrons in 
these insulators. Furthermore, one might note that the 
periodic Anderson model can also be mapped on a Hamil- 
tonian of the form Ti. s + Ti. s f in the limiting case of small 
charge fluctuations. 

Sometimes, especially if higher electron concentrations 
are considered, it is reasonable to include correlations 
within the itinerant electron partial system. The simplest 
possibility is that of an on-site Coulomb repulsion 
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We are interested in the composed Hamiltonian 
"H = TL S + "Hss + tisf + 7~Ljf , 



(5) 



(6) 



which we shall call correlated Kondo lattice model 

(CKLM). By making certain assumptions for the parame- 
ters , U, Jij and J one might return to somewhat more 
special situations. 

The CKLM is not exactly soluble. One needs approxi- 
mation schemes, which are often difficult to justify. Lim- 
iting cases are a powerful tool to weaken this shortcom- 
ing. 
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Till now only a few exact statements on the CKLM 
are known. One of them concerns the limiting case of a 
single excess electron in an otherwise empty conduction 
band, interacting with a ferromagnetically saturated sys- 
tem of localized spins (T = OK). Under these assump- 
tions the election may exist in a bound state, called mag- 
netic polarorMa, which for J < can be shown to be the 
ground stateO. 

Another known and often used exact result for the 
CKLM is that of the zero-bandwidth limitQ. Here one 
assumes the energy dispersion to be flat: e(k) — > T , and 
is then able to calculate all possible one-particle exci- 
tation energies together with the corresponding spectral 
weights. Due to the relation (||) the assumption is equiv- 
alent to Tij — * SijTo, meaning that the a priori itinerant 
electrons described by Hs+Hss are not allowed to change 
the lattice site. For this reason the zero-bandwidth limit 
describes a one-site situation, sometimes called atomic 
limit. 

With this paper we intend to improve the atomic limit 
result by going one step further towards the lattice. A 
two-site cluster possesses two main advantages: First, 
a non-zero intersite hopping integral allows a motion of 
electrons and leads to a finite energy dispersion. Com- 
pared to the flat energy band of the atomic limit, this is 
a qualitative improvement. Secondly, at least two lattice 
sites are necessary to provide a limiting case for an anti- 
ferromagnet. After the two-site CKLM has been solved 
exactly, we suppose to have a powerful tool to test ap- 
proximation schemes for antiferromagnetic solutions of 
the full CKLM (|). Furthermore, we hope that the 
knowledge about properties and symmetries we gained 
when dealing with the two-site problem as well as the 
techniques we used to solve it are helpful tools to develop 
such approximate schemes. 

Especially for the Hubbard model 7i s + 7i ss there have 
already been some attempts to generalize exact cluster- 
results to the lalticcOffllij and to solve even four-site 
clusters exactlyU. In this context one should also men- 
tion efforts to perform a cluster CPA for the Hubbard 
model. 

For the CKLM a complete analytical solution for the 
one-partkis Green's function is so far missing. Only Mat- 
lak et a/£3 have dealt with the eigenvalue problem of a 
two-site cluster for 5 = 1/2 and gave some results for 
certain correlation functions. With our work we tried to 
tackle the task of the exact solution for a two-site cluster 
more comprehensively. The present article is devoted to 
the situation of an insulator, being characterized by an 
empty conduction band. It is our intention to explain the 
employed mathematical techniques and to discuss inten- 
sively the difficulty of a connection between cluster and 
lattice. In a forthcoming articlal3 we will generalize the 
calculations to arbitrary band occupations. 



II. CLUSTER MODEL 

The Hamiltonian (^) possesses with 7i s and 7i// two 
non-local terms. We mentioned before that a two-site 
cluster has two remarkable properties: it is sufficient for 
our intension to improve the atomic-limit result and it 
is still analytically manageable. These desired properties 
guided us to treat the non-local terms as follows: 

The hopping integral is supposed to be 

!T i ^ j and i,j in same cluster, 
To i=j, 
i,j in different clusters, 

leading to a situation as shown in figure Yu. 
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FIG. 1: Formation of two-site clusters due to a restriction of 
the electron hopping. Different symbols for the lattice sites 
are chosen to adumbrate an ABAB-ordered antiferromagnet. 

In our calculations, the Heisenberg term Tiff is ex- 
cluded from the Hamiltonian. This results in a quasi sep- 
aration of the spin system and the electronic system. The 
calculation of the dynamics of localized spins is placed 
back in favour of a description of the conduction band. 
However, a specification of the first influences the latter. 
This offers a chance to connect the two-site cluster to the 
lattice. A localized magnetic moment, calculated within 
the cluster, would give a vanishing value. If (anti-ferro- 
magnetic substances are the subject of investigation, one 
should instead specify the values for (Sf) and (Sf ) to 
be finite. This is only possible if the Hamiltonian of the 
cluster itself does not include a Heisenberg interaction of 
the localized spins. This assumption about our model 
is also justified since in many papers the Kondo-lattiec 
model is understood to consist only oiTL s +'H s f anywayn. 
Matlak et al&3 have chosen another possibility by con- 
sidering a mean-field approximated Heisenberg term as 
part of the Hamiltonian. Then one is confronted with 
even more difficult mathematical expressions than given 
here. 

A simplification we make in this particular article re- 
gards the number of electrons in the conduction band 
N = J2 n i&- As mentioned before, we assume N = 
to be able to give a concise illustration of mathematical 
techniques and physical properties. As a consequence, all 
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expectation values which, like (n^), (^fVr c iT c u)' ' ' 
possess in normal order an annihilation operator are zero. 
However, a treatment with Green's functions implies the 
existence of an additional "test electron" , still giving rise 
to non-trivial results. Nevertheless, the term 7i ss of the 
CKLM, which describes a Coulomb-interaction of two 
conduction electrons, is meaningless in the limit N = 0. 
Therefore, it can be omitted. 

Altogether, our cluster model is described by the op- 
erator ( S 1 * = S+, S l = S~,zi = +1, z l = -1 ): 



Here, we used the Fourier-transformed notation 



n = f E ( c Lc 2ff + 4c lCT ) (?) 

Because it is a finite system, any physical quantity can 
in principle be calculated exactly. Without any restric- 
tion for the concentration of conduction electrons, the 
dimension of the Hilbert space would be 4 2 • (25* + l) 2 . 
We present here results for S = 1/2. For this special case 
one can make use of the operator identity (S?) = h 2 S 2 l, 
which allows to replace any product of two spin operators 
at the same site by a single one: 

st a s? a = h 2 s±hz a s*, (8) 

Sf'S? = -SfSf 17 = TfiSz a Sf a . (9) 

Comparable relations do also existEl for £ > 1/2, giving 
a hint that the qualitative structure of our results is also 
valid for other values of S. 

Even so we are able to solve the eigenvalue problem 
for every electron concentration, it is our intension to 
handle the cluster with a many-particle approach. The 
effort to determine the complete single-particle excita- 
tion spectrum of the cluster is comparable for both tech- 
niques. Additionally, the calculation of the one-electron 
Green's function provides information on the dependence 
of the spectral weights not only on the model parameters 
but also on certain expectation values as the averaged 
z-component of the localized spins. However, the main 
reason for our choice is the wish to get a deeper under- 
standing in the physics and possible approximations of 
the full CKLM, where a many-particle treatment is in- 
evitable. 



III. SOLUTION FOR THE EQUATIONS OF 
MOTION 

Due to the Heisenberg time dependence of the opera- 
tors the retarded one-particle Green's function 



(10) 



A(t);B(t'))) :=-ie(t-t')^[A(t),B(t') 
has the following equation of motion 



E((A;B)) E = h([A,B]. 



[a,hub)) e . (ii) 



A:B 



d(t-t') ((A(t);B(t'))) exp 



h 



E(t - t') 



Motivated by the site symmetry, it is convenient to 
consider the following combinations of Green's functions: 



CsD- s c sa 



(<0 



G s D s c S!J ] c sa 



G s D s c S(J \ c sa 



(12) 



Here, the site index s represents the opposite of site s 
(s = 1 =>■ s = 2; s = 2 s = 1), C s and D s are arbitrary 
products of spin operators at site s and s, respectively, 
e = +1 or —1. 

We are especially interested in the one-electron 



Green's function 



The com- 



plete set of Green's functions generated by a repeated 



imp 



emcntation of (h if) is 
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G ( 7 e) 
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Kh 


= o, 




= «sr^O) (e) 


K 9 


= n(sr°sz) 




= ((S s a SgC sa )) 


K w 





(13) 

In the right column the corresponding inhomogeneities 
in the equation of motion ( |ll] ) to each of the Green's 
functions is given. Emerging Green's functions with an- 
other combination of spin operators or a higher number 
of electron operators vanish identically because of the 
constraints S — 1/2 and N = 0. 



It is a valuable property of the combination (12) that 
the equations of motion of the G^ and the G^ Green's 
functions do not mix. Moreover, if one writes ( p"l| ) in the 
form 
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z J m ij Gj — Ki , 



(14) 



one gets a matrix M — (n r iij) i ■ which contains e only as 

a pre- factor of T. If the matrix equation (fi4) is solved 

for e = +1, one obtains the expression for G\ simply by 

replacing in g\ + ^ the parameter T by —f. We therefore 
omit e in the forthcoming equations. 
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By a couple of unitary transformations M can be re- 
duced further to blocks of a size not greater than 3x3. 
It is favourable to introduce the combined Green's func- 
tions 



H, 



(±) 



h 

2 \\^S<7 



((c s J±z a ((S z s c s J 



H 



(±) 



I {(S§c s J ± z a ((S* s S§c s J 



H. 



(±) 



lfr(±) 



2 ^rr-'-'i 



(±) 



^3 - 2^0 ±^1 



= I {(St"S?"S§c a J 



(± 



4 - 2 ((^X-.)) ± z a ({SIS-°c ±s _ J 

(±) _ 



(15) 

In another step the combinations ' ± G10, #3 ± Gg 
and H[ +) ± are used, what finishes the reduction. 

The remaining 3x3 blocks have the following matrix- 
structure: 



M, 



(/Jl/X2A'3A 1 4) 



(I 



E + 2/i4a fiiT 
(E) = I Ml T £ ^ 2 2a 
/i 2 2a E + fj, 3 f 



Here, we used the abbreviation a 



(16) 



2 2 and ^1, ... ,^4 
are sign parameters. The eigenvalues of this matrix are 



(M3AU) 

02 

(/J3/J4) 

03 



4a 2 - 2n 3 fi 4 af + f 2 , 
Aa? - 2 f j, 3 ii i af + f 2 . 



For the eigenvectors one gets 



(^1^2/^3/^4) 
01 



1 / MlM3 
' 1 



and for k = 2, 3 



v, 



(mi/^MsAu) 
Ofe 



V3 



Pk 



^1 (T - 2// 3 ^ 4 a - fi^Eok) 
2/i 2 a 



with 



1 E ak +n 3 T-fi 4 a 



2a 



3-Eofc 



(17) 



(18) 



(19) 



We denote by VQk[i] the ith component of Vq^ and de- 
fine a matrix V by its components Vo[ifc] = «ofe[i]. Ba- 
sic algebra shows that Mq = Vq ■ Dq ■ Vq with Dq = 
diag(i? + E ok ). Using these expressions one can give the 
components of the inverse of Mq in a form which is char- 
acterized by linear energy poles: 



3 

M o l [ij] =^2v ok [i]v k\j] 



1 



E + E k 



(20) 



This form is particularly suitable for the calculation of 
Green's functions. 

The 10 x 10 matrix equation (|lj) can now be written 
in the form 



M, 



(-+— 



M, 



(-+-+) 



(E-T Q 





Mt +++) (E-T -a) 



and solved simply by applying equation (|2(j) . The emerg- 
ing inhomogeneities R\ , . . . , Rg are constructed out of 
K±, . . . , -ftTio, obeying the same rules of combination as 
for the Green's functions. The single missing equation 
has an even simpler structure. One can write down im- 



(+). 



mediately the solution for H 2 

H ( 2 +) = (21) 

E-T a -T + a 

§ (1ft 2 + Tiz a (Sf)) + (\tfz a (Sj) + h (SjSfi) 
E-To-f + a 

To obtain an expression for the one-particle Green's 
function we are looking for, one has to take a sum of the 



combined functions 



,if 5 (±) . Hence, 



- Gi - -o (H. 



r(+) 



H. 



(+) 



H. 



(-) 



H. 



(-) 



(22) 



— ^ I H, 



r(+) 



H. 



(-) 



+ 
+ 



2?r 



H. 



H. 



(-) 
3 

(-) 



G10 
G10) + ( 



H. 



H. 



(+) 
3 

(+) 



Gg 

Gg 



We now remind the fact stated after equation (jlj) that 
this solution can be obtained equally well for e = +1 as 
for e = — 1. If we add these two Green's functions 

«c SCT ;cL» E = \ («0> (+) + <(^)) ( " } ) , (23) 

we effectively reverse the combination (^), which was a 
consequence of site symmetry. 

Therefore, the following result for the one-particle 
Green's function is obtained: 
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£-£3 
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E -e 4 



E-£6 



un 2 



(-2z a h(S z ) + z a h{S z y ) + 2(S Z S Z ) 



h h h 



E-e 3 
1 



E - 



E-e 4 



E 



h 2 - ZrT h(s z x ) ~ z a h(s z y ) 
-2(S Z S Z ) (SfS*) - (s-«s%)) 



h/&- +) h/E^ +) h/E { ++y 



^03 



E-e 3 
f 



E — £5 E — £4 E — £q 



It has the six energy poles 



(24) 




Site 1 



0.25 0.5 0.5 0.25 

Hybridization T [eV] 



Site 2 



FIG. 2: The hybridization T causes a splitting of the one- 
electron energy levels valid for a single site and a group- 
ing of the excitation energies. For T = 0.5 eV and Jh = 
0.3eV,Tb = OeV the arrangement of the energy poles E = 
£2, £5, £6, £i, £4, £3 (starting from the bottom) is shown. 
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A. Remaining in the cluster 



As a result of the calculations in section III we do 
not only get an expression for the Green's function G± 
and therefore for ((c scr ; cJ CT ))„. On the contrary, we are 
without further effort able to write down similar expres- 
sions for all the other participating Green's functions 



G 



functions H, 



One just has to sum the combined Green's 

(±) w(±) 

5 



, Hk^' another way as in (C 



It is one possible concept to apply the well-known spec- 
tral theorem to each of the determined Green's functions. 



If G[ e) , 



G 



(<0 



or G±Q is of the form 



IV. TREATMENT OF THE CORRELATION 
FUNCTIONS 



A-B 



E+io+ 



£ 



E - Ei + i0+ ' 



(25) 



The result (^4J) provides an expression for the one- 
particle Green's function. It is characterized by six 
linear energy poles in the denominators. They are 
identical with the excitation energies possible if a sin- 
gle excess electron is placed in the otherwise empty 
conduction band. The numerators give the spec- 
tral weights corresponding to these excitations. How- 
ever, in contrast to the model parameters, T, Tq,J, 
which are assumed to be known, we have till now 
no information about the spin correlation functions 
(Sf) , (Si) , (SfS$) , (S+S 2 -) , (S^S+), leaving the spec- 
tral weights undetermined. 

In the next two subsections we will introduce and dis- 
cuss two possible concepts to obtain results for these cor- 
relation functions. They are especially important for the 
treatment of a situation with N > 0. 



then one can derive an expression for a corresponding 
correlation function 

(si) = = £ e(3iE Z + 1 ' ( 26 ) 

i=i i=i 

with the help of this theorem (/3 = (fee? 1 ) -1 ). One gets 
a set of ten equations, which are linear in the unknown 
correlation functions. To solve this systems means that 
only properties of the cluster determine the correlations 
functions. 

However, in our special situation of N — each of the 
expectation values at the left hand side of (^6|) is sup- 
posed to vanish, because they all include an electron den- 
sity operator. Then the solution of the set of equations 
becomes trivial. One has to choose \i <C £fcV7c, to fulfill 
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the requirements for the electron density. This leaves the 
set of spin correlation functions again undefined. 

Hence, if JV = 0, the application of the spectral 
theorem to the determined (that means electronic) set 
of Green's functions fails to provide any information 
about spin correlation functions. Instead other Green's 
functions can be used. The equations of motion for 
((S+; Sf» B , <(5+; Sj)) E and «5|5+; Sf» B are partic- 
ularly simple in the case N = 0. If the spectral theorem 
is applied to its results, one obtains the information that 
(Sf ) , (Sf) and (SfS%) are zero, respectively. However, 
even when using ^Sf; S^~)) B or ((S^~; Sf)) E it is not 

possible to get any expression for (S^S^) and (S^S^)- 
For these correlation functions the expectation value has 
to be determined directly by calculating the trace of the 
canonical ensemble with the four possible eigenstates of 
H for N = 0. The result is again zero. For an insulator it 
is not surprising that all spin correlation functions vanish 
if they are solely determined within the cluster. 



B. Connection to the lattice 

To create a connection of the two-site cluster to the 
lattice, the spin correlation functions need to be specified 
outside the cluster. It is the advantage of an insulator 
that this can be done without coming into conflict with 
the spectral theorem applied to the cluster. 

As explained in section [ij] we have quasi separated the 
electronic system and the spin system by excluding Tiff 
from the Hamiltonian. As a different problem, we can 
now consider a situation of two localized spins mutually 
interacting via TL ff with = Jh and embedded in the 
lattice via an additional term in the Hamiltonian of the 
form 



Hf = -b(SZ + riS$) 



(27) 



Simultaneous to the choice of Jh, the sign parameter 
r\ allows to distinguish between ferromagnetic (rj = +1) 
and antiferromagnet (77 = —1) configurations. In the first 
case b — B ■ 9J £ B is as a molecular or Weiss field, in the 
latter case it describes the crystal anisotropy. The com- 
plete set of Green's functions for this proble m is de rived 
in appendix JaI As a result the expressions (A13)-(A16) 
and ( A21 )-( A24 ) for a set of correlation functions are 
achieved. 

This set is identical with the correlation functions we 
need for the calculation of the electronic part described 
by TL. Determined in this way and used in (B4), we are 



sure that not only (Sf) and (S|) but also the two-site 
correlation functions (Sf Sf) , (S^ Sf") and (SrS^~) are 
described consistently by two experimentally accessible 
parameters Jh and B. We created an intelligible connec- 
tion to the lattice. 

Only after these considerations we are able to specify 
the behaviour of the electrons in a two-site (anti-ferro- 
magnetic Kondo lattice. If we assume a ferromagnet 



{(SI) — (Sf)), we obtain a dependence of the spectral 
weights belonging to the six energy poles of the single- 
electron Green's function on the magnetization as shown 
in figure [| The magnetic field B is uniquely used to cre- 




Ma^iLUK' field B [T] 



FIG. 3: For the energy poles of the single-electron Green's 
function, the dependence of the corresponding spectral 
weights on a magnetization (Sf) = (Sf) parallel to the spin 
of the test electron a is given (left). The latter is adjusted 
by an assumed molecular field B (right). The fixed model 
parameters are: J H h 2 = 10" 3 eV, Jh = 0.3 eV, f = 0.5 eV 
and g.j = l,To = OeV. The temperature is 80 K. 

ate a ferromagnetic order of the system of localized spins 
and does not affect the conduction electrons directly. To 
obtain a saturation of about 95% one needs extremely 
large exchange fields of about B = 400 T, as can also be 
estimated if the tendency to disorder, k^T ', and to order, 
HbB, are set equal. One can clearly see that for these 
fields the energy poles £2 and £1 become dominant. This 
is not surprising, since the configuration where the ex- 
cess electron and the two localized spins have all parallel 
spin has exactly this excitation energy. The observation 
that for B = T these two poles still possess the highest 
spectral weight is connected to the fact that addition- 
ally some one-electron eigenstates of the Hamiltonian TL 
with the z-component of the total spin ±(2S — 1/2) have 
the same energy. The symmetry of lattice sites 1 and 2 
nicely enters the figure, too: It is the reason for similar 
structures in the groups around E = —T = — 0.5eV and 
around E = f = 0.5 eV. 

As one can see in figure || the situation is quite dif- 
ferent if the spin of the test electron is antiparallel to 
the alignment of the localized spins, induced by the ap- 
plied magnetic field. Here, the energy poles £3 and e§ 
get maximum spectral weight when the field is increased. 
However, the site symmetry is not affected. 

If an antiferromagnetic configuration of the localized 
spins ((Sf) = — (Sf)) is considered, there is again no 
change in the positions of the energy poles (figure |^). 
Of course this does not apply to their spectral weights. 
Comparing figure || to figure [| one notices already for 
zero magnetization ((Sf) = <^ B = 0) a different dis- 
tribution. The opposite sign for Jh affects the non-zero 
correlation functions (SfS%) and (Sf^S^)- As concerns 
the dependence on B, primarily the behaviour of the en- 



FIG. 4: Same as figure g with the only difference that the 
spin of the test electron a is antiparallel to the direction of 
the magnetization. 



FIG. 6: Same as figure |5| but the Green's function is evalu- 
ated at the site which has a localized spin antiparallel to the 
electron spin. 




FIG. 5: Same as figure §, but with (Si) = - (Sf) and J H h 2 = 
-l(T 3 eV. B is the magnetic anisotropy. 



energy levels of the two single sites, as shown in figure |2j. 
The fact that the spectral weights of the energies in the 
lower group get the values of the higher group (and vice 
versa) if the spin direction of the test electron is changed, 
is again related to site symmetry. 

Of course, all these features can be explained by a de- 
tailed analysis of the eigenvalue problem. However, the 
Hilbert spaces that need to be considered are relatively 
complex. For instance the subspace corresponding to a 
z-component of the total spin of 1/2 is 6-dimensional, 
and hence every eigenvector is a linear combination of 
6 different spin configurations. Therefore, the possible 
analysis is too lengthy to be presented here. 



ergy poles e$ and £4 has changed. Their spectral weight 
does not tend to zero (ferromagnet), but rather increases 
with B. In general the antiferromagnetic configuration 
leads to smaller changes of the spectral weights, since the 
average \ ((Sf) + (<Sf )) = remains constant. 

The term TLf 1 which enters the determination of the 
correlation functions, leads for r\ = —1 to a break-down 
of the site symmetry. Therefore, the two groups of en- 
ergy poles in figure a show a slightly different dependence 



E 

The 



on B. One has to distinguish between ucj-pC^ 

z n> c ti)) ? and (( c w c n)) E = (( c 2f> c h)) E 

latter situation has been presented in figure 0. The de- 
pendence on the spin direction of the electron is of course 
much smaller than for a ferromagnetic configuration. In 
fact, apart from an interchange of the two groups of en- 
ergy poles the distribution of the spectral weights shown 
in figure || is almost the same as in figure |[ The forma- 
tion of these two groups is a result of hybridization of the 



V. SUMMARY 

We have shown that the Hamiltonian Ti can be treated 
such that an analytic expression for the one-particle 
Green's function as given in ( p4| ) is obtained. In order to 
do that we formulated the complete set of equations of 
motion as a single matrix equation, made use of the sym- 
metries of the problem and combined Green's functions 
to reduce the matrix to blocks of size 3x3, which can be 
solved. In this article we studied the case of an insulator 
(N = 0), but a generalization of these techniques to the 
situation in a metal is possible even so calculations are 
much less concise there. i3 

Furthermore, the connection to a lattice has been dis- 
cussed. We propose to use the correlation functions in the 
spectral weights of the Green's functions for this purpose. 
A helpful tool is a two-site Heisenberg model which has 
been separated before from the Hamiltonian (^|) and cal- 
culated in appendix |A]. By adding a term Tif we can sim- 
ulate ferromagnetic and antiferromagnetic lattices. This 
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assumption about the system of localized spins results 
in certain features of the electronic partial system, visu- 
alized exemplary by the spectral weights given in figure 
^ to ^. The site symmetry has a major influence on its 
distribution. 
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APPENDIX A: THE TWO-SPIN PROBLEM 

It is the purpose of this appendix to provide a consistent set of spin correlation functions as input for the two-site 
CKLM. This is done by discussing the Heisenberg-type Hamiltonian 

n A = - j h (s+s^ + s^s+ + 2Sj£|) - b (si + v si) . (Ai) 

To obtain correlation functions we use the spectral theorem, which needs certain Green's functions as an input. Within 



an equation-of-motion (11) approach these are 



G n = ((S^;S+)) E K n = -2h 2 (S[), 
Gai = ((S 2 ;S+)) E K 21 = 0, 

r 12 = ((s* 1 Sz;S+)) E l 12 = n 2 (s+s 2 ), 



Following the rules identical to what has been used in ([lfj) the simplification 

E + bh±j H n 2 T-m 2 \ /|G n ±r 21 \ /|K n ±i 21 

TJnh 2 E + ijbh± j H n 2 J \ §g 2 i ± r 12 J \±k 21 ±l 12 

is possible. The determinant of the emerging matrix, has the form 

1 1 / 1 1 



(A2) 



T 21 = «5|5r;5+)) B L 21 = -2ft 2 (SfSf/ • 

On the right hand site we give the corresponding inhomogeneities which appear in the equations of motion: 

[E + bK)G n = K 11 + 2J H h{T 12 -r 21 ) (A3) 
{E + v bh)G 21 = K 21 -2J H h(T 12 -r 21 ) (A4) 

h 3 

(E + r)bK)T u = L 12 + J H y (G n - G 21 ) (A5) 

h 3 

(E + bfi)Y 21 = L 21 - J H y (Gn - G 2i ) (A6) 



(A7) 



D(±) 2y^T ys± J H n 2 + 1(1 + 77)6-^ E± J H ?i 2 + |(l+7?)6 + 
with J7~. = J(J u h 2 ) 2 + ^(l-7 1 )b^h 2 . (A8) 



By taking the inverse of the matrix in (A7) and after performing a partial fraction expansion one therefore obtains 
the combined Green's functions 

2 Gn±r 2 i — -h (S r 1 ) 7-t 7-7 — 

V E-Ef' E-E { 2 ±} J 
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T 2^ 2 SfSf ^ 



( S+s -) Jh ^ 4 [ 1 1 
|G 21 ± r 12 = ± (5 ?> ^ ( - ^-^j 



J H ft 4 / 1 1 



2 / 



2 



E-E^ E-E^ 



with the four energy poles 



= T J H a 2 -|(l + r 7 )6--y/(j H /i 2 ) 2 + ^(l-?7)^ 2 , 



4 ±} = Trf- |(1 + r,)6+^/(j H ?i 2 ) 2 + i(l- t?)62^. (A9) 
One can apply the spectral theorem ( p6|) to these Green's functions. The corresponding left hand sides are 

G 2 i ► (Si~S 2 )) C Ain ^ 

r 12 — (5+5f5 2 -) = -1(5+^), (A10) 

r 2 i — (Sfs+sr) = ^(si) + h(sfsi). 



The emerging Bose distribution functions at the right hand side of ([260 can ^ e abbreviated: 



™f 3 = ^ • ( A11 ) 



c 



per - l 



One ends up with a system of three equations which are linear in the correlation functions we are looking for. The 
fourth equation is just necessary to calculate It is most instructive to give the corresponding matrix equations 

for each of the following two special situations separately. 



1. Ferromagnet (77 = +1) 

Here, the expressions for the energy poles are 



(+) 
1 

(+) _ 



-hb~2J H h 2 , e{~ 



= -fib, 



E. 



= -fib, 

) = -hb + 2J H h 2 . 



(A12) 



As a consequence m\ ' = . The system of equations is: 




-1 



±m 3 



— if m\ + ^ — mi, 



I (mi 



+) 



,(-) 



5 777,4 



where the abbreviations 777,3 — ntJi + 2m\ 1 + vn\ ' and 7774 = 777^ 



,(-) 



,(-) 



4 +) 



-1 



,(+) 



m\ — 777 



2 ^"'1 

1 + im 3 



(-) 



,(+) 



2777 



(-) 



,(-) 



h(Sf) 

2 (Sf Si) 
(StSz) 



have been used. 
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This [gives a solution for the correlation functions which can be simplified to the expressions published before by 



Nolti 



<s+s 2 -) 

(SIS Z 2 ) 
(SI) 



exp[/3K,&] — exp[— (3hb] 



2 1 + exp[-2/?fr 2 J H ] + e^q>[f3hb] - exp[-/37i&] ' 



1 - exp[-2/m 2 J H ] 



2 1 + exp[-2/3?i 2 JH] + exp\j3hb] - exp[-/?7i&]' 
n exp [(3hb\ - exp[-f3hb] - 1 - exp[-2/3?i 2 JH] 

4 1 + exp[-2/3^ 2 J H ] + exp[/?n&] - exp[-0tO>]' 
(SI) 



In contrast to his calculations we obtain here the last line as a result and not as an assumption. 



(A13) 

(A14) 

(A15) 
(A16) 



2. Antiferromagnet (rj = — 1) 

Here, the expressions for the energy poles are 



e{ +) = 


-./ H ?i 2 - 




K^) 2 > 


(A17) 


£f> = 


+ J H ?i 2 - 






(A18) 


4 +) = 


- J H ?i 2 + 


\/(rf) 2 - 




(A19) 


4-> = 


+ J H ?i 2 + 


V/W 2 ) 2 - 


h (bhf 


(A20) 



and consequently 



4 +) 
,(+) 



mi, } = -1, 

m{ ' = — 1, 



(+) (-) o (+) 
m ( 2 +) - m { f ] = 2m 2 +) 



1, 
1. 



The system of equations is: 



o I = 





,(+) 



,(+) 



,(+) 



(+) 



,(+) 



,(+) 



4 +) ) 



V 



'2 



7= H 



+) 



l 2 



2 (Sf 5|) 



7 / 



Here, the following set correlation functions, which to our knowledge has not been published before, is obtained. 



{st s^) = 



(S{S Z 2 



(St) = 



,(+) 



,(+) 



(jW) +(bh) 



1 2 + Am^m^ + 



m[ +) + m 2 +) 



bh 



m 



(+) 



rn 



(+) 



y/ '(Jnh 2 ) 2 + (bh) 

(si) = - (st) 



2 2 + Am\ + ' > m 2 +S> + m 2 + ^ 



(A21) 

(A22) 

(A23) 
(A24) 



The last line, which characterizes an antiferromagnet, is again a result of the calculations. 
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